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ANALYSIS  OF  THE  HESSIAN  FOR  INVERSE  SCATTERING 
PROBLEMS.  PART  III:  INVERSE  MEDIUM  SCATTERING  OF 
ELECTROMAGNETIC  WAVES  IN  THREE  DIMENSIONS 


TAN  BUI-THANH  *  AND  OMAR  GHATTAS  t 

Abstract.  Continuing  our  previous  work  [6,  Inverse  Problems,  2012,  28,  055002]  and  [5,  In¬ 
verse  Problems,  2012,  28,  055001],  we  address  the  ill-posedness  of  the  inverse  scattering  problem  of 
electromagnetic  waves  due  to  an  inhomogeneous  medium  by  studying  the  Hessian  of  the  data  misfit. 
We  derive  and  analyze  the  Hessian  in  both  Holder  and  Sobolev  spaces.  Using  an  integral  equation 
approach  based  on  Newton  potential  theory  and  compact  embeddings  in  Holder  and  Sobolev  spaces, 
we  show  that  the  Hessian  can  be  decomposed  into  three  components,  all  of  which  are  shown  to 
be  compact  operators.  The  implication  of  the  compactness  of  the  Hessian  is  that  for  small  data 
noise  and  model  error,  the  discrete  Hessian  can  be  approximated  by  a  low-rank  matrix.  This  in 
turn  enables  fast  solution  of  an  appropriately  regularized  inverse  problem,  as  well  as  Gaussian-based 
quantification  of  uncertainty  in  the  estimated  inhomogeneity. 

Key  words.  Inverse  medium  scattering;  electromagnetic  wave  propagation,  adjoint,  Hessian, 
Newton  potential,  compact  operators,  compact  embeddings,  Riesz- Fredholm  theory,  ill-posedness, 
Gauss-Newton. 

AMS  subject  classifications.  49N45,  49K20,  49K40,  49K30,  49J50,  49N60 


1.  Introduction.  A  feature  of  many  ill-posed  inverse  problems  is  that  the  Hes¬ 
sian  operator  of  the  data  misfit  functional  is  a  compact  operator  with  rapidly  decaying 
eigenvalues.  This  is  a  manifestation  of  the  typically  sparse  observations,  which  are 
informative  about  a  limited  number  of  modes  of  the  infinite  dimensional  field  we  seek 
to  infer.  The  Hessian  operator  (and  its  finite  dimensional  discretization)  play  an  im¬ 
portant  role  in  the  analysis  and  solution  of  the  inverse  problem.  In  particular,  the 
spectrum  of  the  Hessian  at  the  solution  of  the  inverse  problem  determines  the  degree  of 
ill-posedness  and  provides  intuition  on  the  construction  of  appropriate  regularization 
strategies.  This  has  been  observed,  analyzed,  and  exploited  in  several  applications 
including  shape  optimization  [12,13]  and  inverse  wave  propagation  [4,11,14],  to  name 
a  few. 

Moreover,  solution  of  the  inverse  problem  by  the  gold  standard  iterative  method — 
Newton’s  method — requires  “inversion”  of  the  Hessian  at  each  iteration.  Compactness 
of  the  Hessian  of  the  data  misfit  functional  accompanied  by  sufficiently  fast  eigenvalue 
decay  permits  a  low  rank  approximation,  which  in  turn  facilitates  rapid  inversion 
or  preconditioning  of  the  regularized  Hessian  [8,11].  Alternatively,  solution  of  the 
linear  system  arising  at  each  Newton  iteration  by  a  conjugate  gradient  method  can 
be  very  fast  if  the  data  misfit  Hessian  is  compact  with  rapidly  decaying  eigenvalues 
and  the  conjugate  gradient  iteration  is  preconditioned  by  the  regularization  operator 
[1],  Finally,  under  a  Gaussian  approximation  to  the  Bayesian  solution  of  the  inverse 
problem,  the  covariance  of  the  posterior  probability  distribution  is  given  by  the  inverse 
of  the  Hessian  of  the  negative  log  likelihood  function.  For  Gaussian  data  noise  and 
model  error,  this  Hessian  is  given  by  an  appropriately  weighted  Hessian  of  the  data 
misfit  operator,  e.g.,  [20].  Here  again,  exploiting  the  low-rank  character  of  the  data 
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misfit  component  of  the  Hessian  is  critical  for  rapidly  approximating  its  inverse,  and 
hence  the  uncertainty  in  the  inverse  solution  [4,7,14,17]. 

In  all  of  the  cases  described  above,  compactness  of  the  data  misfit  Hessian  is  a 
critical  feature  that  enables  fast  solution  of  the  inverse  problem,  scalability  of  solvers  to 
high  dimensions,  and  estimation  of  uncertainty  in  the  solution.  With  this  motivation, 
here  we  analyze  the  Hessian  operator  for  inverse  medium  electromagnetic  scattering 
problems,  and  study  its  compactness.  This  is  a  continuation  of  our  previous  efforts 
in  analyzing  the  ill-posedness,  via  the  compactness  of  the  Hessian,  of  inverse  shape 
acoustic  scattering  [5]  and  inverse  medium  acoustic  scattering  [6]  problems.  Our 
analysis  is  based  on  an  integral  equation  formulation  of  the  time  harmonic  Maxwell 
equations,  adjoint  methods,  and  compact  embeddings  in  Holder  and  Sobolev  spaces. 
These  tools  allow  us  to  analyze  the  Hessian  of  the  data  misfit  in  detail.  The  ideas 
of  proof  are  similar  to  those  in  [6],  but  the  details  are  more  technical  due  to  the 
complexity  of  the  Maxwell  equations.  Numerical  results  demonstrating  the  theoretical 
results  are  not  carried  out  in  this  paper  since  similar  numerical  validations  can  be 
found  in  our  previous  work  [5,6]. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2  briefly  derives  and 
formulates  forward  and  inverse  electromagnetic  scattering  problems  due  to  bounded 
inhomogeneity.  We  then  derive  the  Hessian  for  the  inverse  problem  in  Section  3. 
Section  4  justifies  the  Hessian  derivation  by  studying  the  well-posedness  of  the  (in¬ 
cremental)  forward  and  (incremental)  adjoint  equations,  and  the  regularity  of  their 
solutions  with  respect  to  the  smoothness  of  the  inhomogeneity.  Next,  we  analyze  the 
Hessian  in  Holder  spaces  in  Section  5,  and  then  extend  the  analysis  to  Sobolev  spaces 
in  Section  6.  Finally,  the  conclusions  of  the  paper  are  presented  in  Section  7. 

2.  Forward  and  inverse  medium  problems  for  electromagnetic  scatter¬ 
ing.  In  this  section,  we  briefly  discuss  forward  acoustic  scattering  problems  due  to 
bounded  inhomogeneity  and  the  corresponding  inverse  problems.  We  begin  by  assum¬ 
ing  that  the  time  harmonic  incident  wave  (E1C,  HJC)  satisfies  the  Maxwell  equation 

VxEic  —  ikHic  =  0,  inR3, 

VxH“  +  ik~Eic  =  0,  inR3, 

where  k  >  0  is  the  wave  number,  i2  =  —  1  the  imaginary  unit.  The  total  field 
formulation  for  the  scattering  due  to  dielectric  medium  [10]  can  be  equivalently  cast 
into  the  following  scattered  field  formulation: 

VxE-  jfcH  =  0,  in  R3,  (2.1a) 

VxH  +  iknE  =  ik(l  —  n)Eic,  in  R3,  (2.1b) 

lim  (H  x  x  —  rE)  =  0,  r  =  ||x||  ,  (2.1c) 

r—¥  oo 

where  n  =  n  (x)  >  0  is  the  distributed  refractive  index  which  is  assumed  to  be  1  for 
the  free  space  (if  not,  we  can  always  normalize  (2.1a)-(2.1b)  to  fulfill  this  assumption), 
and  (E,  H)  the  scattered  field.  The  Silver-Miiller  radiation  condition  (2.1c)  is  assumed 
to  be  valid  uniformly  in  all  directions  with  x  G  R3  denoting  the  vector  of  spatial 
coordinates.  Note  that  it  can  also  be  written  as  limr.  r  (H  x  n  —  E)  =  0,  where 
n  is  the  unit  outward  normal  vector  of  the  sphere  with  radius  r.  To  the  rest  of  the 
paper,  the  medium  is  assumed  to  be  bounded,  that  is,  there  exists  a  bounded  domain 

H  such  that  n(x)  =  1,  V  ||x||  eR3\Sl.  In  other  words,  q  =  1  —  n  has  compact  support 
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in  fi.  In  fact,  we  shall  assume  n  £  C'm’“(R3),  and  hence  q  £  throughout 

the  paper,  where  m  £  N. 

For  the  forward  problem,  n  is  given  and  we  solve  the  forward  equations  (2.1a)- 
(2.1c)  for  the  scattered  field  (E,  H).  For  the  inverse  problem,  on  the  other  hand,  given 
observation  data  (Eobs,Hobs)  over  some  compact  subset  flobs  C  M3,  we  are  asked  to 
infer  the  distribution  of  the  refractive  index  n.  One  way  to  solve  the  inverse  problem 
is  to  cast  it  into  the  following  PDE-constrained  optimization  problem: 

min  J  :=  [  K{x)  f  |E  -  Eobs|2  +  |H  -  Hobs|2N)  dx,  (2.2) 

«  Jr  a  V  / 

subject  to  the  forward  equations  (2.1a)-(2.1c).  Here,  K(x)  is  the  observation  operator 
whose  support  is  Qobs.  Note  that  we  could  have  different  observation  operators  for  E 
and  H  but,  for  simplicity  of  the  exposition,  we  assume  they  are  the  same  and  equal 
to  K(x).  In  order  to  cover  several  interesting  observation  operators,  fiobs  is  allowed 
to  be  quite  general  in  this  paper.  In  particular,  it  could  be  a  closed  subset  in  M3  or  a 
relative  closed  subset  of  a  manifold  in  M3.  For  example,  Hobs  could  be  a  closed  arc, 
or  a  closed  curved,  or  a  closed  subset  of  a  two  dimensional  manifold,  or  some  two 
dimensional  manifold.  For  convenience,  we  identify 

ICip:=  /  Kip  dx  =  /  ip(y)dy. 

J  R3  J  n°b‘ 

We  also  permit  pointwise  observations  in  our  analysis,  i.e.  flobs  =  ,  and  in 

this  case  we  identify 


Kip  := 


(2.3) 


3.  Derivation  of  the  Hessian.  In  this  section,  we  derive  the  gradient  and 

Hessian  of  the  data  misfit  (2.2)  using  a  reduced  space  approach1,  and  the  justification 

for  our  derivations  is  provided  in  Section  4.  We  begin  with  a  useful  observation  on 

the  radiation  condition.  Since  the  radiation  condition  (2.1c)  is  valid  uniformly  in  all 

directions  A,  we  rewrite  the  radiation  condition  as 
11*11  ’ 


H  x  n  —  E  =  ip(r)  =  o  (r  3)  , 

where  r  is  the  radius  of  a  sufficiently  large  circle  Too,  and  ip(r)  =  o  (r^1)  means 


lim  rip(r)  =  0. 


It  can  be  seen  that  the  cost  functional  (2.2)  is  real-valued  while  the  constraints 
(2.1a)-(2.1c)  are  complex- valued.  Consequently,  the  usual  Lagrangian  approach  will 
not  make  sense  and  care  must  be  taken.  Following  Kreutz-Delgado  [16],  we  define  the 


1See  [19]  for  the  difference  between  the  full  space  and  the  reduced  space  approaches. 
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Lagrangian  as 

C  =  J~j  h  •  ( V  x  E  —  ikH)  dx 

J  R3 

+  J  er  ■  [H  x  n  —  E  —  ip(r)]  ds 
-  [  h-  (VxE  +  ifcH)  dx  +  [ 


+  /  e  •  [VxH  +  iknE  -  ik(  1  -  n) Eic]  dx 

J  R3 


/  R3 


VxH  —  iknE  +  ik{  1  —  ?^)E^ 


dx 


+ 


er  ■  [H  x  n  —  E  —  y?(r)]  ds 


where  the  overline,  when  acting  on  forward  and  adjoint  states  (and  their  variations), 
denotes  the  complex  conjugate. 

Taking  the  first  variation  of  the  Lagrangian  with  respect  to  e,  h,  er  in  the  direc¬ 
tions  e,h,er  and  arguing  that  the  variations  e.h,  er  are  arbitrary  yield  the  forward 
equations  (2.1a)-  (2.1c). 

Now  taking  the  first  variation  of  the  Lagrangian  with  respect  to  E,  H  in  the 
directions  E,  H,  using  the  following  vector  calculus  identities 


V  •  (A  x  B)  =  B  ■  (Vx A)  -  A  •  (VxB) , 

A  •  (B  x  C)  =  B  ■  (C  x  A)  =  C  •  (A  x  B) , 


(3.1a) 

(3.1b) 


and  arguing  that  the  variations  E,  H  are  arbitrary  yield  the  following  adjoint  equa¬ 
tions: 


and 


Vxe  —  ikh  =  —K  (H  -  Hobs)  , 
Vxh  +  ikne  =  K  (E  -  Eobs)  , 
lim  (h  x  x  +  re)  =  0, 

r— too 

er  =  —  h  x  n  on  ] 


m  K 
in  R: 
r  =  llxll 


(3.2a) 

(3.2b) 

(3.2c) 


It  should  be  pointed  out  that  the  adjoint  equations  (3.2a)-(3.2c)  is  very  similar  to  the 
forward  equations  (2.1a)-(2.1c)  except  for  the  minus  sign  in  the  radiation  condition. 
This  can  be  understood  that  the  adjoint  waves  originate  from  infinity. 

If  we  eliminate  er,  the  Lagrangian  now  becomes 

C  =  J~l  h-  (VxE-AH)  dx+  /  e  •  [V x H  +  iknE  —  ik(l  —  n)E*c]  dx 


/R3 


/R3 


(h  x  n)  •  [H  x  n  —  E  —  <p(r)]  ds 


-  h-  (VxE  +  ifcH)  dx 


/R3 


VxH  —  iknE  +  ik(  1  —  n)Etc 


dx 


J  (h  x  n)  •  [H  x  n  —  E  —  <p(r)]  ds 


The  gradient  of  the  cost  function  acting  in  the  direction  q  is  simply  the  variation  of 
the  Lagrangian  with  respect  to  q  in  the  direction  q ,  i.e. , 


DJ  (q;q)  =  ik 


/  R3 


(e  +  EIC) 


-e  -  ( E  +  E*' 


gdx, 


(3.3) 
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which  is  clearly  a  real  number. 

For  the  sake  of  convenience  in  deriving  the  Hessian,  the  forward  and  adjoint 
equations  are  best  expressed  in  the  weak  form.  As  a  direct  consequence  of  the  above 
variational  calculus  steps,  the  forward  equation  in  the  weak  form  reads 


S  (q,  E,  H)  =  -  /  h  •  (VxE  —  ikU)  dx  +  /  e  •  [VxH  +  ikrtE  -  ik{  1  -  n)EiC]  dx 


—  j  ^xnj'[Hxn-E  -  ¥5(r)]  ds  =  0,  V  (e,  hj  . 
Similarly,  the  adjoint  equation  in  the  weak  form  is  given  by 

A  (q,  E,  H,  e,  h)  =  [  H  •  [Vxe  -  ikh  +  K  (H  -  Hobs)]  dx 

J  R3 

-  [  E  •  [Vxh  +  ikne  -  I\  (E  -  Eofcs)]  dx 


(3.4) 


/  R3 


x  •  (h  x  n 


e)  ds  =  0,  V  (  E,  H 


(e,h) 


(3.5) 


Next,  the  (reduced)  Hessian  acting  in  the  directions  q  and  q  is  obtained  by  simply 
taking  the  first  variation  of  the  gradient  DJ  ( q ,  q)  with  respect  to  q ,  E,  H,  e  and  h  in 
the  directions  q,  E,  H,  e  and  h,  i.e. , 

D2J{q-q,q)=ik  [  [e  •  (e  +  E1^)  -  d  •  (E  +  Eic)  +  e  •  E  -  e  •  e]  q  dx  (3.6) 
J R3  L  v  /  J 

As  mentioned  at  the  beginning  of  this  section,  the  reduced  space  approach  is  employed, 
and  hence  the  variations  E  and  e  can  not  be  arbitrary.  In  fact,  they  are  only  admissible 
if  the  forward  and  adjoint  equations  are  satisfied.  As  a  direct  consequence,  the  first 


variations  of  S  (q,  E,  H)  and  A  (q,  E,  H,  e,  h)  must  vanish,  that  is,  ^E, 
solution  of  the  following  incremental  forward  equation: 


is  the 


/  11 '  ( 


VxE  —  ikHj  dx  + 
(hxn).(Hxii-E) 


VxH  +  iknE  —  ikq  (E  +  E* 


dx 


ds  =  0,  V 


(e,h). 


(3.7) 


and 


(e,  hj  i 

J  5-( 


is  the  solution  of  the  following  incremental  adjoint  equation: 


Vxe  —  ikh  +  /AH  )  dx  — 


/  E'( 


V  x  h  +  ikne  —  K E  —  ikqe )  dx 


J  (h  x  n)  •  (h  x  n  +  e)  ds  =  0,  V  (e, h)  . 


(3.8) 


Consequently,  the  corresponding  strong  form  of  the  incremental  forward  equation  is 

in  R3,  (3.9a) 

in  R3,  (3.9b) 


VxE  —  ikH.  =  0, 

VxH  +  ikriE  =  ikq  (E  +  Ei! 
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and  that  of  the  incremental  adjoint  equation  reads 

Vxe  —  ikh  =  -A'H, 

Vxh  +  ikne  =  A'E  +  ikqe , 

lim  hxx  +  re  =0, 
r— >oo  \  J 


in  R3,  (3.10a) 

in  R3,  (3.10b) 

r  =  ||x||.  (3.10c) 


It  should  be  pointed  out  that  the  incremental  forward  pair  ^E,  HJ  in  (3.9a)-(3.9c)  is 
a  function  of  q,  i.e. ,  ^E  (g) ,  H  (g)^j  for  a  given  variation  q.  Similar  observation  applies 
to  the  incremental  adjoint  pair  (e(g),h(g)J  in  (3.10a)-(3.10c).  To  avoid  possible 


confusion,  if  E  is  evaluated  at  q,  for  example,  we  will  write  this  dependency  explicitly 
as  E  (q) ;  otherwise  E  means  E  (q)  implicitly. 

Next,  we  need  to  convert  the  Hessian  in  (3.6)  into  a  symmetric  form  that  is 
convenient  for  our  later  analysis.  The  first  step  is  to  replace  q  by  q  and  choose 
^e,  =  ^e(g),h(g)^  in  the  incremental  forward  equation  (3.7).  In  the  second  step, 

we  take  ^E,  =  ^E(g),H(g)^  in  the  incremental  adjoint  equation  (3.8).  The 

last  step  is  to  subtract  the  resulting  incremental  forward  equation  from  the  complex 
conjugate  of  the  resulting  incremental  adjoint  equation.  After  some  simple  integration 
by  parts  using  (3.1a)  and  cancellations,  we  obtain 


ik 


e  (q)  ■  (e  +  Eic')  —  e  (q)  ■  (E  +  EIC)  qdx  =  ik  e  ■  E  (q)  —  e  ■  E  (q)  qdx 

V  /  ']  ./ TB3  L 

I<  [H  (q)  ■  H  (q)  +  H  (q)  ■  H  (5)1  dx  +  [  K 


E  (q)  ■  E  (q)  +  E  (q)  ■  E  (q) 


dx 


which,  after  combining  with  (3.6),  gives  the  desired  symmetric  form  of  the  Hessian  as 


D2J(q;q,q)=  /  K  E  (q)  ■  E  (q)  +  E  (q)  ■  E  (q) 


dx 


I< 


Wi  (q\q,q) 

U(q)-W{q)+W{q)-U(q) 


dx 


'Hiki'AA) 


ik 


)  R3  L 


e  ■  E  (q)  -  e  •  E  (q) 


qdx  +  ik 


e  ■  E  (q)  -  e  •  E  (q) 


(3.11) 


qdx . 


4.  Regularity  of  the  forward  and  adjoint  solutions.  In  this  section  we 
are  going  to  justify  what  we  have  done  in  Section  3  by  studying  the  regularity  of 
the  (incremental)  forward  and  (incremental)  adjoint  solutions  with  respect  to  the 
medium  q  and  its  variations  g,  q.  In  particular,  for  sufficiently  smooth  inhomogeneity, 
the  solutions  turn  out  to  be  classical  by  using  an  integral  equation  method,  as  we 
shall  show.  As  a  by-product,  we  shall  also  show  that  the  (incremental)  forward  and 
(incremental)  adjoint  equations  are  well-posed. 

Let  us  introduce  the  following  standard  volume  potentials  (also  known  as  Newton 
potentials)  [10,15]: 


w(x)=Ttp(x)=  [  $(x,y)v?(y)dy,  x  G  R3, 

Jr3 


(4.1) 
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where  4>  is  either  the  fundamental  solution  of  the  (incremental)  forward  equation(s) 
defined  as 


4>(x,y) 


eife||x-yJ[l 

4?r  ||x  —  y ||  ’ 


or  the  fundamental  solution  of  the  (incremental)  adjoint  solution(s): 


$(x,y) 


e— ifc||x— y|| 

4tt  ||x  —  y || ' 


We  recall  the  standard  mapping  properties  of  T  defined  in  (4.1)  here. 

Lemma  4.1.  Let  ip  £  Co  (R3)  D  (^""’“(R3),  where  a  £  (0, 1]  and  m  £  N.  Then 
w  £  Cm+2’a  (fi),  and  ||w||Cm+2,Q(n)  <  c\\ip\\Cm,a^,  where  supp  (ip)  C  O  C  R3. 
Moreover,  the  orders  of  differentiation  and  integration  can  he  interchanged. 

Proof.  See  [6,  Lemma  1]  for  a  proof.  □ 

The  following  representation  formula  is  an  easy  variant  of  the  Stratton  and  Chu 
theorem  [10,  Theorem  6.1]. 

Lemma  4.2.  Let  £  C 1  (R3),  there  holds 


£  (x)  =  -  lim  V  x  /  n  (y)  x  S  (y)  (x,  y)  ds  (y) 

r->oo  JTao 

+  lim  v/  n(y)  •  (y)  <J>  (x,  y)  ds(y) 

r-s-oo  JTao 

-  lim  ik  n(y)x  Jf  (y)  (x,  y)  ds  (y) 

r->oo  JT ^ 

+  Vx  f  [Vx^  (y)  —  ikJ£  (y)]  4>(x,  y)  dy 

J  R3 


-V  /  V  •  #(y)<L(x,y)dy 

J  R3 


+  ik  [Vx^(y)  +  ik£  (y)]  $(x,  y)  dy, 

J  R3 

Vx  e  R3. 

(4.2) 

Next,  we  show  that  the  boundary  terms  in  (4.2)  vanish  for  a 

large  class  of 

problems. 

Lemma  4.3.  Let  £,Jff  £  C 1  (R3)  satisfy 

\7x£  -ikJ£  =  f, 

in  R3  \  fl, 

(4.3a) 

V  xJF  +  ik£  =  g, 

in  R3  \  11, 

(4.3b) 

lim  (J^  x  x  ±  r£)  =  0, 

r—too 

r  =  ||x|| , 

(4.3c) 

23?  f  Jtf  •  f  ds  <  oo,  23?  f  £  ■  gds  <  oo, 

(4.3d) 

Jr „ 


where  3?  is  an  operator  that  returns  the  real  part  of  its  argument.  Then 

lim  [  \\Jif  x  n|| 2  ds  =  O  (1) ,  and  lim  f  \\£\\2  ds  =  O  (1) . 

r—foo  J p  r—foo  J p 
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Moreover, 


-  lim  Vx  [  n  (y)  x  £  (y)  $  (x,  y)  ds  (y) 

r->oo  7r=o 

+  lim  V  /  n  (y)  ■  £  (y)  $  (x,  y)  ds  (y) 

r-s-oo 

—  lim  ik  n  (y)  x  Jif  (y)  <!>  (x,  y)  ds  (y)  =0,  Vx  £  l3. 

r->oo  J  roo 

Proof.  First,  from  the  assumption  on  the  radiation  condition  (4.3c)  we  have 

J  \\j£  x  n±#||2  ds  =  0, 


from  which  it  follows  that 

x  n|jz  ± 


ds  =  q=23?  /  (n  x  £)  ■  ds. 


Second,  integrating  both  sides  of  (3.1a)  on  K3  \  with  A  =  #,  B  =  Jif  and  using 
Gauss  divergence  theorem  together  with  assumptions  (4.3a)-(4.3b)  give 


/  (n  x  £)  ■  ds+  (n  x  £)  ■  ds  =  ik 
J  Jan  J  r 


/R3\fi 


/R3\n 


dx 


£  ■  g  dx  ±  /  ■  f  dx. 

J  H3\n 


It  follows  that 


=F23?  /  (n  x  £)-j£ds  =  ±23?  /  (n  x  £)-Jif  ds±2iR  /  #-gdx=F23?  /  Jif-fdx. 


/an 


l  R3\n 


Consequently,  we  have  the  identity 

\\j£  x  n|| 2  +  ||# ||2)  ds  =  ±23?  /  (n  x£)-j£ds 


ion 
±23? 


/  £  ■  g  dx  =f  23?  /  •  f  dx, 

/R3\n  dR3\n 


whose  right  side  is  finite,  and  so  is  the  left  side.  The  first  assertion  is  therefore 
justified.  The  second  assertion  is  the  direct  consequence  of  the  first  and  the  proof 
of  [10,  Theorem  6.6].  □ 

Define,  for  any  xell, 

«(y) 


(Q[q\£)  (x)  =  V 


fn  1  -  g(y) 


V  •  #(y)$(x,y)dy, 


K(£,J£)(x)  = -Vx  I  A'(y)  jT4>(x,y)dy  ±  •  f  - — 1--d\  (y)#$(x,  y)  dy 

*  dn^  1  -  g(y) 

±*fc  f  AT(y)#$(x,y)dy, 

jQobs 

(T[q}£)  (x)  =  V  /  - - !--Vg  (y)  ■  #  (y)  $(x,  y)  dy  ±  k2  f  q( y)#(y)$(x,  y)  dy, 

dn  i  —  <?  fyj  Jn 

(T[g]#)  (x)  =  V  /  - - ^77r:yVd(y)  -#(y)$(x,y)dy±fc2  [  q(y)£(y)$(x,  y)  dy, 


i-«(y) 
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and  I  as  the  identity  operator.  Here  is  a  direct  consequence  of  Lemma  4.1. 

Corollary  4.4.  Suppose  q,  q  £  C™’a(£l),  m  £  N. 

i)  Q  maps  Cs+1’^  (f2)  continuously  into  Cs+1't3  (0)  for  s  +  (3  <  m  +  a. 

ii)  Cartesian  components  oflZ  are  analytic. 

in)  T,T  map  (f2)  continuously  into  Cp,a  (fl)  for  1  <  p  <  to,  and  T,‘T  : 

Cp~  1,“  (fi)  — >  CS'P  (f2)  are  compact  operators  for  p— l+a<s  +  )3<p  +  a. 

Proof.  It  is  sufficient  to  prove  the  third  assertion.  A  straightforward  application 
of  Lemma  4.1  shows  that  T  is  a  continuous  operator  from  Cp~1,a  (fi)  to  Cp’a  (fi). 
Moreover,  the  compactness  is  trivial  due  to  compact  embeddings  in  Holder  spaces  [21]. 
□ 

We  are  in  the  position  to  discuss  representation  formulas  for  (incremental)  forward 
and  (incremental)  adjoint  equations. 

Proposition  4.5.  Solutions  of  the  (incremental)  forward  and  (incremental) 
adjoint  equations  satisfy: 


(T  +  T[q])  E  =  —T[q]  Eic, 

(4.4) 

(J  +  T  [q])  E  =  -f  [q]  (E  +  Elc)  -  Q[q]  E, 

(4.5) 

(l  +  T[q])e  =  K(  E  -  Eobs,  H  -  Hobs)  , 

(4.6) 

(1  +  T  [?])  e  =  n  (E,  H)  -  t(q\e  -  Q[g]e. 

(4.7) 

Conversely,  suppose  that  E,E,  e,  e  are  solutions  of  the  integral  equations  (4.4)-(4.7), 
and  define,  Vx  €  K3, 


E(x)  :=  -T  [q]E-T  [q]  Eic, 

E(x)  :=  -T  [q]E-T  [g]  (E  +  EIC)  -  QE, 
e(x)  :=  -T  [q\  e  +  K  (E  -  Eobs,  H  -  Hobs)  , 
e(x)  :=  -T  [q]  e  +  K  (E,  h)  -  f  [q}e  -  Qe, 

H  :=  VxE/(ifc),  h  :=  Vxe/(*fc)  +  K  (H  -  Hobs)  /(ik),  H  :=  VxE/(ifc),  and  h  := 
Vxe/(ife)  +  KH/(ik).  Then  (E,  H) ,  (e,h)  ,  (e,h)  ,  ^e,  are  solutions  of  the  (in¬ 
cremental)  forward  and  (incremental)  adjoint  equations. 

Proof.  We  first  observe  that  the  (incremental)  forward  and  (incremental)  adjoint 
equations  satisfy  Lemma  4.3  since  q,  q  and  I\  have  compact  support.  Lemma  4.2  then 
shows  that  all  these  equations  have  the  following  common  representation  formula 

#(x)=Vx  [  [V x#  (y)  —  ikJi?  (y)]  <I>(x,  y)  dy 

J  R3 

-V  /  V-<T(y)$(x,y)dy 

J  R3 

+  ik  I  [VxJ^(y)  +  ik$(yj\  <I>(x,  y)dy,  Vx  e  l3.  (4.8) 

Jr3 

We  provide  the  proof  for  only  (4.7)  since  the  others  are  simpler.  From  (3.10b)  we  find 
that 

V  ■  e  =  —  —  Vn  •  e  — •  rj*-V  •  (a"e)  +  —V  •  (qe) , 
n  kn  \  /  n 
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which,  together  with  (3.10a)-(3.10b),  can  be  substituted  to  (4.8)  to  end  the  proof  of 
(4.7)  by  using  the  fact 


/  -V 

/  R3  n 


(A'E)4>(x,y)ciy 


K(y ) 

n(y) 


E(y)$(x,y)dy, 


owing  to  Vx4>(x,  y)  =  —  Vy$(x,  y),  suppqnsuppiv  =  0,  and  integration  by  parts. 
The  converse  is  long  and  tedious,  and  hence  omitted,  though  one  can  similarly  follow 
[10,  Theorem  9.2].  □ 

Proposition  4.5  shows  that  we  can  equivalently  study  the  existence,  uniqueness, 
and  stability  of  (incremental)  forward  and  (incremental)  adjoint  solutions  via  the 
integral  equations  (4.4)-(4.7). 

Theorem  4.6.  Let  q,q  G  (7™’a(12),  m  £N,  there  hold: 
i)  there  exists  a  unique  solution  to  the  forward  integral  equation  (4.4)  such  that 


e e  cm’a(n),  He  cm+1 

II H 


I  cm+u 


qn)  -  c  I 


(O) 


{fl),  <md||E||c. 
"E  +  E’ 


qn)  —  c  I 


’(«) 


E*' 


ii)  there  exists  a  unique  solution  to  the  incremental  forward  integral  equation 


(4.5)  such  that  EEC" 


where 


E 


C" 


<  C 


(12),  H  G  Cm+1'a  (12),  and 

IE 


H 


C™  +  V 


\Cn 


qn) 


I C” 


Ei( 


I C" 


qn) 


y 


<  c 


E 


in)  there  exists  a  unique  solution  to  the  adjoint  integral  equation  (4.6)  such  that 


e  G  Cm,a  (f2)  and 


IC" 


qo) 


<  c  (||E  -  E 


obs  I 


H  H 


obs  I 


ii)  there  exists  a  unique  solution  to  the  incremental  adjoint  integral  equation  (4.7) 


such  that  e  G  Cm’a  (11)  and 


qn)  -  c  1 


E 


H  +  llellcm,“(n) 

oo  y 

is  analytic  Corollary  4.4 


Proof.  We  begin  with  the  first  assertion.  Since  E1 
shows  that  the  right  side  of  (4.4)  belongs  to  Cm’a  (12),  and  that  T  is  compact  in 
Qrn-1,0  gy  the  Riesz-Fredholm  theory  [9],  we  only  need  to  prove  the  uniqueness 
of  (4.4).  That  is,  we  need  to  prove  that  (I  +  T[q\)  E  =  0  in  C™_1,a  (12)  implies 
E  =  0  in  Cm~1,a  (12).  Now,  Proposition  4.5  indicates  that  E  and  H  =  VxE/(ifc)  is 
the  solution  of  (2.1a)-(2.1c)  with  E*c  =  0.  Following  Step  2  of  the  proof  of  Lemma 
4.3  in  12  with  S  =  E,  ^  =  H,  f  =  g  =  0,  we  find  that 


2 3?  /  (n  x  S’)  ■  Jf?ds  =  0, 

Jon 

which,  by  Steps  1  and  2  of  the  proof  of  Lemma  4.3,  implies  that 

lim  [  ||E||2  ds  =  0. 

r-J-ooJr^ 

Now,  since  (E,  H)  satisfies  the  Maxwell  equation  in  R3  \  12,  each  component  of  E 
satisfies  the  Helmholtz  equation  in  R3  \  12,  and  we  conclude  that  E  =  H  =  OinR3\12 
by  Rellich  lemma  [10].  A  result  in  [10,  Theorem  9.3]  then  shows  that  E  =  H  =  0  in  R3. 
Invoking  the  Riesz-Fredholm  theory  yields  that  there  exists  a  unique  E  G  Cm~1’a  (12) 
and 


El 


qn) 


<  cIlTME" 


qn) 


<  c| 


qn) 


E*' 


Icm-ncqn) 


where  we  have  used  [5,  Lemma  5]  in  the  last  inequality.  Unless  otherwise  stated, 
the  constant  c  at  different  places  may  have  different  value.  Using  Lemma  4.1  shows 
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that  in  fact  E  £  (7m’a  (f2)  and  ||E||CmiC,(a^  <  c||E||Cm_i,„^.  Next,  using  Lemma 
4.1  and  the  fact  VxV  (•)  =  0,  we  see  that  H  £  Cm+1’a  (LI)  and  < 

cll9llc"*'“(fi)  llE  + 

Similarly,  for  the  incremental  forward  equation,  we  observe  that  the  right  side 
of  (4.5)  lives  in  Cm,a  (f 2)  by  an  application  of  Lemma  4.1.  An  application  of  the 
Riesz-Fredholm  theory  shows  that  E  £  Cm~ 1,a  with 

|E||0„.,,.(n)  <  «ll«llc-..(0)  (l|E|lc...,n)  +  ||E"||c.,.(n))  , 

and  the  property  of  the  Newton  potential  in  Lemma  4.1  improves  the  result  to  E  £ 
Cm’a  {n). 

The  proofs  for  the  adjoint  and  incremental  adjoint  equations  follow  similarly  by 
observing  that  1Z  (•,  ■)  is  analytic  due  to  the  fact  Llobs  n  fl  =  0. 

□ 

Remark  4.7.  Let  us  define  min  {(m,  a) ,  (s  +  1,  /3) }  to  be  equal  to  (m,  a)  if  m  + 
a  <  S  +  1+/3  and  to  (s  +  l,/3)  otherwise.  Then,  a  more  general  result  can  be  obtained: 
for  the  forward  equation,  for  example,  if  1  <  p  <  m,  for  p  —  l  +  a<s  +  /3<p  +  a 
there  exists  a  unique  solution  to  (4.4)  such  that  E  £  Cmin {(jn,a),(s+i,p)}  jj  g 

Cmin{(m+i,a),(s+2^)}(0)_  Moreover,  ||E||Cmin{(m,Q)i(5+1^}(n)  <  c  ||Eic||Cm_1|Q(n),  and 

1 1 n ||gmiii{(m+i,c,),(s+2,/3)}  <  c  ||E*c||gm_i,c,(f2)-  However,  for  clarity  of  the  exposition,  we 
consider  only  the  case  p  =  s  =  m ,  /3  =  a. 

We  are  now  in  the  position  to  justify  our  derivations  of  gradient  and  Hessian  in 
Section  3. 

Theorem  4.8.  Let  Q  C  R3  be  a  bounded  domain.  Assume  that  q,  q  and  q 
belong  to  Cg’a(f2).  Then,  the  cost  functional  (2.2)  is  twice  continuously  Frechet 
differentiable;  and  hence  the  gradient  (3.3)  and  Hessian  (3.11)  are  well-defined. 

Proof.  First,  observe  that  we  have  used  Gateaux  derivatives  to  derive  the  gradient 
and  Hessian  in  Section  3.  Now  it  is  evident  that  both  DJ(q;  q)  and  D2J (q;  q ,  q)  are 
linear  and  continuous  with  respect  to  q  (and  q)  by  Theorem  4.6.  Moreover,  continuous 
dependence  on  q  of  E,  H  from  Theorem  4.6  implies  the  continuous  dependence  on  q  of 
e,  E,  and  H,  which  in  turn  implies  the  continuity  of  DJ(q\q)  and  D2J(q]q,q)  with 
respect  to  q.  Hence,  a  classical  result  on  sufficiency  for  Frechet  derivative  [3]  ends  the 
proof.  □ 

5.  Analysis  of  the  Hessian  in  Holder  spaces.  In  this  section  we  study  the 
behavior  of  the  Hessian  at  a  fixed  refractive  distribution  n,  i.e. ,  q  =  1  —  n  £  C(rfl'°  (fl). 
Similar  to  our  previous  work  [6]  on  medium  inverse  acoustic  scattering,  the  Hessian 
of  the  inverse  medium  scattering  problem  for  electromagnetic  waves  turns  out  to  be 
compact  for  all  q  as  we  shall  show.  For  concreteness,  we  restricted  ourselves  to  two 
exemplary  cases  of  the  observation  operator,  namely,  the  observation  is  everywhere  on 
a  compact  subset  having  non-trivial  r-dimensional  Lebesgue  measure  for  some 
1  <  r  <  3  (we  call  this  case  as  continuous  observation)  and  pointwise  observation 

nobs  =  {xfs}^\ 

From  Theorem  4.6,  observe  that  the  incremental  forward  solution  E  can  be  iden¬ 
tified  as  the  following  operator  composition: 

C0m’Q  (fi)  3  n  ^  E (q)  =  -  (1  +  T[(?])_1  T[q]  (E  +  EIC)-(I  +  T[g])_1  QE  £  Cm~1,a(Ll). 
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Lemma  5.1.  E  (g)  :  C™’Q  (ft)  — >  Cm  1,a(f2)  is  a  compact  operator. 

Proof.  First,  owing  to  the  Riesz-Fredholm  theory,  (/  +  T[g])_1  is  a  bounded 
map  from  to  Cm-^a(ft).  Second,  due  to  E  £  Cm’a  (ft)  C  Cm-1<a(ft), 

Corollary  4.4  shows  that  both  T  and  Q,  as  a  function  of  g,  map  (ft)  continuously 
into  Cm,a  (ft),  and  hence  compactly  embedded  into  Cm~1,a(ft)  [21] .  Finally,  since 
E  is  a  sum  of  two  operators,  each  of  which  is  a  composition  of  a  continuous  and  a 
compact  operators,  it  is  compact.  □ 


As  a  direct  consequence,  E(-) 


n°b» 


:  (ft)  ->•  C m-1’Q  ( flobs )  is  still  a  compact 


operator  since  restricting  to  flobs  is  a  continuous  operation. 

If  the  observation  is  continuous,  the  first  Gauss-Newton  part  of  the  Hessian, 
namely  TLi(q~,  q ,  q),  can  be  now  rewritten  as 


H1(q;q,q)  =  2M  (E(n),E(n)) 


L2(Oofes) 


=  23?  (  E 


E(n) 


,h)  ,  , 

nobs  /(C“'“(fi))xc0"*'°(n) 


where  (•)*  denotes  the  adjoint  operator.  Note  that  the  composition  E 


E 


IS 


meaningful  since  Cm  1,a  ( flobs )  C  ( Cm  1,01  (ftobs))';  in  fact  it  is  continuous.  In  this 
form,  the  first  part  of  the  Gauss-Newton  Hessian 


7*1  (?)  (■:•) 


H1(g;-,-)  =  2R 


* 


n°b- 


.  To  study  the  second  part 


is  evidently  compact  due  to  the  compactness  of  E  (•) 
of  the  Gauss-Newton  Hessian,  we  first  write  H  (q)  as 

H  (q)=  ^-VxE(q)  =  ikf7x  [  g(y)E(g)$(x,y)dy+ifcVx  f  q(y)  (E  +  E*c)  $(x,  y)  dy, 
ik,  Jo  Jo 


where  we  have  used  the  fact  that  VxV(-)  =  0.  Since  E  (q)  maps  (H)  con¬ 
tinuously  into  Cm,a  (f2),  we  conclude  that  H  (q),  as  a  operator,  is  continuous  from 
C™’0  (ft)  to  Cm+1’a  (H)  owing  to  E  £  Cm,a  (H),  q,  q  £  C™’Q  (ft),  and  Lemma  4.1.  By 
the  standard  Holder  compact  embeddings  [21],  H  (q)  is  a  compact  operator  from  from 
C™’a  (ft)  to  from  Cm,a  (ft).  Argue  similarly  as  above,  we  deduce  that  the  second  part 
of  the  Gauss-Newton  Hessian 


%(g)(v)  :=^2(g;-,-)  =  25R(fi  H(-) 


Qobs 


is  a  compact  operator. 

Let  us  now  study  the  Gauss-Newton  Hessian  with  pointwise  observations.  On  the 
one  hand,  we  rewrite  the  right  side  of  (4.5)  as 


-  V, 


7y.(g(E  +  Eic)) 


(E  +  E1' 


$(x,y) 
n(  y) 
$*(x.y) 


dy-k2  [  g  (E  +  Eic)  $(x,  y)  dy 
J  n 


»(y) 


—  k2  (E  +  Elc)  $(x,y) 


q(  y)  dy, 


fi(x,y) 


where  integrating  the  first  term  by  parts,  using  compact  support  of  q,  and  interchang¬ 
ing  the  order  of  integration  and  differentiation  have  been  used.  Here,  =  Vx^  and 
Vy  is  understood  to  act  component-wise. 
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On  the  other  hand,  using  the  adjoint  E*  allows  us  to  write 


T[q]  E  (x)  =  -  (  E* 


^(x,y) 
i  -  9  (y) 


v?  (y)  ,  q( y)  )  -  (  E*  {k2q( y)$(x,  y)) ,  q{ y) 


f2(x,y) 


b(x,y) 


with  (-,•)  denoting  the  duality  pairing  between  (C™’Q  {Vi))'  and  C™'a  (12).  Conse¬ 
quently,  from  incremental  forward  integral  equation  (4.5),  the  evaluation  of  E  (q)  at 
x*  can  be  written  as 


E(g)  (xf )  =  ([fr(xf  ,y)  -  f2(xf  ,y)  -  f3(xf s, y)] , g(y)^  . 

It  follows  that  the  first  part  of  the  Gauss-Newton  Hessian  can  be  written  as 

Nobs  _  /  /  Nobs  \  \ 

Hi(q)  (q,  q)  =  23?  £  E  (q)  (xf >)  E  (g)  (xf >)  =  ^  E  ^  ,  9^  , 

which  shows  that  the  dimension  of  the  range  of  T-Li(q)  is  at  most  Nobs.  As  a  result, 
7~li(q)  is  a  compact  operator.  By  the  same  token,  one  can  show  that  the  second  part  of 
the  Gauss-Newton  Hessian  T~Li(q)  (<?,  q)  is  compact  due  to  its  finite  dimensional  range. 

We  summarize  the  above  result  on  the  compactness  of  the  Gauss-Newton  Hessian 
in  the  following  theorem  which  is  valid  for  both  continuous  and  pointwise  observation 
cases. 

Theorem  5.2.  The  Gauss-Newton  Hessian,  TLi(q)  plus  TL2 (q),  as  a  continuous 
bilinear  form  on  Cff  ,a  (0)  x  Cq1’“  (f2),  is  compact. 

The  analysis  of  Tl^q)  is  somewhat  easier  as  we  shall  now  show. 

Theorem  5.3.  TL3  (q),  as  a  continuous  bilinear  form  on  C™’“  (f l)  x  C™’a  (fl),  is 
a  compact  operator. 

Proof.  Rewrite  7^3(9;  q,  q)  as 


in  L 


'H3(q;q,q)  =  2kiR  ie  ■  E  (q)q  +  ie  ■  E  (q)q 


dft  =  2 m  (ie  ■  E  {q)  +  E*  (ieq)  ,  q 


We  conclude  that  TL3  (q)  is  compact  by  the  following  three  observations.  First,  the 
incremental  forward  solution  E  can  be  identified  as  a  compact  operator  as  discussed 
above.  Second,  multiplication  by  ie  is  a  continuous  operation  in  Cm’a  (Q)  (see  [5]  for 
example).  Third,  the  sum  of  two  compact  operators  is  again  compact.  □ 

We  close  this  section  by  observing  that  the  full  Hessian  is  the  sum  of  three  compact 
operators,  it  is  therefore  compact  as  well. 


6.  Analysis  of  the  Hessian  in  Sobolev  spaces.  Similar  to  our  previous  work 
[5,6],  we  shall  extend  the  analysis  in  Holder  spaces  to  Sobolev  spaces.  We  recall  a 
result  [6]  on  the  mapping  property  of  T  in  Sobolev  spaces  that  is  result  similar  to  that 
of  Lemma  4.1. 

Lemma  6.1.  Assume  that  is  bounded  and  integrable,  fl  C  M3  is  a  bounded 
domain.  Then  T  defined  in  (4.1)  maps  Hm  (12)  continuously  to  Hm+2  (12)  for  m  £ 
NU{0}. 

We  next  summarize  main  results  parallel  to  those  in  Sections  4  and  5.  Since  the 
proofs  are  similar,  we  omit  the  details. 

Corollary  6.2.  Suppose  9, 9  £  Hm  (12),  to  £  N. 
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i)  Q  maps  Hm  (12)  continuously  into  Hm  (12). 

ii)  T,  T  map  iLm_1  (12)  continuously  into  Hm  ( 12),  and  T,  T  :  Hm~1  ( 12)  — > 
Hp  (12)  are  compact  operators  for  m  —  1  <  p  <  m. 

Note  that  the  last  result  in  Corollary  6.2  is  due  to  compact  embeddings  in  Sobolev 
spaces  [2, 18].  From  now  on  to  the  end  of  this  section,  we  need  to  >  2  for  point 
observations  to  make  sense  since  only  Sobolev  spaces  of  order  greater  than  3/2  are 
embedded  into  the  space  of  continuous  functions  [18] . 

Theorem  6.3.  Let  q,q  £  Hm  (12),  TO  £  N  and  to  >  2,  there  hold: 
i)  there  exists  a  unique  solution  to  the  forward  integral  equation  (4.4)  such  that 


E  £ 

IIH 


Hm  (12),  H  £ 

<  cIMI 


Hm+1(  SI) 


Hm+1  (12),  and 

E  +  E 


E 


H 1 


- 


E4' 


Hm(  n) 


h™(  n)' 


ii)  there  exists  a  unique  solution  to  the  incremental  forward  integral  equation 


(4.5)  such  that  E  £  Hm  (12),  H  £  Hm+1  (12),  and 


H 


H™+1(Q) 


<  C 


E 


where 


E 


h™(  n) 


<  c||g|| 


J?m(f 2) 


El 


H™(n) 


+  E4' 


)■ 


Hi)  there  exists  a  unique  solution  to  the  adjoint  integral  equation  (4.6)  such  that 


e  £  Hm  (12)  and 


<  c(||E  -  E( 


obs  I 


H  H 


obs  I 


ii)  there  exists  a  unique  solution  to  the  incremental  adjoint  integral  equation  (4.7) 


'(fi)  —  c  I 


E 


H 


+  e 


H™(Q) 


such  that  e  £  Hm  (12)  and  ||e||ff„ 

Clearly,  the  Sobolev  setting  also  shows  that  the  gradient  and  the  Hessian  are 
meaningful. 

Theorem  6.4.  Let  12  C  R3  be  a  bounded  domain.  Assume  that  q,  q  and  q 
belong  to  Hm  (12)  with  to  £  N  and  m  >  2.  Then,  the  cost  functional  (2.2)  is  twice 
continuously  Frechet  differentiable;  and  hence  the  gradient  (3.3)  and  Hessian  (3.11) 
are  well-defined. 

Lemma  6.5.  Let  q,  q  £  Hm  (12),  to  £  N  and  m  >  2.  Then  E  (g)  :  Hm  (12)  — > 
Hm~l  (12)  is  a  compact  operator. 

Here  comes  the  main  result  of  this  section. 

Theorem  6.6.  Let  12  c  K3  be  a  bounded  domain.  Suppose  q  £  Hm  (12),  to  £  N 
and  m  >  2.  Then,  the  Hessian,  TL  ( q )  =  Hi  ( q)+H2  (<?)+%  (q),  is  a  compact  operator 
in  Hm  (12). 


7.  Conclusions.  We  have  analyzed  the  Hessian  stemming  from  the  inverse  prob¬ 
lem  of  scattering  of  electromagnetic  waves  due  to  bounded  inhomogeneity.  Similar  to 
our  companion  paper  on  inverse  medium  scattering  problems  of  acoustic  waves  [6] ,  we 
have  shown  that  the  full  Hessian  is  a  compact  operator.  Our  analysis  starts  with  a 
study  on  the  regularity  of  the  scattering  solution  with  respect  to  the  smoothness  of  the 
medium  based  on  the  Newton  potential  theory  and  the  Riesz-Fredholm  framework. 
Then,  together  with  compact  embeddings  in  Holder  and  Sobolev  spaces,  we  are  able  to 
prove  the  compactness  of  the  Hessian  operator  in  both  Holder  space  and  Sobolev  space 
settings  for  three  dimensional  inverse  medium  electromagnetic  scattering  problems. 
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